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Abstract 

If future neutrino oscillation experiments show that the neutrino mass spectrum 
is with normal ordering, mi < m2 < m^, and the searches for neutrinoless dou- 
ble beta {{(3P)ou-) decay with sensitivity to values of the effective Majorana mass 
|<m>| 10-2 eV give negative results, the next frontier in the quest for (/3/3)oi/-decay 
will correspond to |<m>| ~ 10"^ eV. Assuming that massive neutrinos are Majorana 
particles and their exchange is the dominant mechanism generating (/?/3)o;y-decay, we 
analise the conditions under which |<m>|, in the case of three neutrino mixing and 
neutrino mass spectrum with normal ordering, would satisfy |<m>| > 0.001 eV. We 
consider the specific cases of i) normal hierarchical neutrino mass spectrum, ii) of rel- 
atively small value of the CHOOZ angle ^13 as well as iii) the general case of spectrum 
with normal ordering, partial hierarchy and a value of ^13 close to the existing upper 
limit. We study the ranges of the lightest neutrino mass mi and/or of sin^^is, for 
which |<m>| > 0.001 eV and discuss the phenomenological implications of such sce- 
narios. We provide also an estimate of |<m>| when the three neutrino masses and 
the neutrino mixing originate from neutrino mass term of Majorana type for the (left- 
handed) flavour neutrinos and '^'j mjU^j = 0, but there does not exist a symmetry 
which forbids the (/3/3)oiy-decay. 
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1 Introduction 



The experiments with solar [H [2], [3], atmospheric [3], reactor [5l E] and accelerator 
neutrinos [7[ |8] have provided during the last several years compelling evidence for the 
existence of neutrino oscillations caused by nonzero neutrino masses and neutrino mixing. 
The neutrino oscillation data (see also [9l |10]) imply the presence of 3- neutrino mixing in 
the weak charged lepton current (see, e.g. [TT]): 



where i^/l are the flavour neutrino fields, Uji^ is the field of neutrino Uj having a mass rrij and 
U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [12], U = f/pMNS- 

In spite of the remarkable progress made, first, in demonstrating experimentally the 
existence of neutrino oscillations and, second, in determining the pattern of neutrino mixing 
and the values of the two neutrino mass squared differences, responsible for the solar and 
atmospheric neutrino oscillations, our knowledge in what concerns most of the basic aspects 
of neutrino mixing is very limited at present (see, e.g. [H]). We still do not know i) what 
the nature of neutrinos with definite mass is - Dirac or Majorana, ii) what type of spectrum 
neutrino masses obey, iii) what the absolute scale of neutrino masses is, iv) whether the 
CP-symmetry is violated in the lepton sector by the neutrino mixing matrix ?7pMNS5 v) what 
the value of the CHOOZ angle is - being the smallest mixing angle in the PMNS matrix, it 
controls (together with the Dirac CP-violating phase) the magnitude of CP-violation effects 
in neutrino oscillations, vi) whether the observed patterns of neutrino mixing is related to 
the existence of a new symmetry in Nature, etc. 

Establishing whether the neutrinos with definite mass uj are Dirac fermions possessing 
distinct antiparticles, or Majorana fermions, i.e. spin 1/2 particles that are identical with 
their antiparticles, is of fundamental importance for making progress in our understanding of 
the origin of neutrino masses and mixing and of the symmetries governing the lepton sector 
of particle interactions (see, e.g. [H]). It is well-known that the presence of massive Dirac 
neutrinos is associated with the existence of a conserved additive lepton number, which can 
be, e.g. the total lepton charge L = Le + + L^-. If the particle interactions do not conserve 
any lepton charge, the massive neutrinos uj will be Majorana fermions (see, e.g. [13]). 

The only feasible experiments having the potential of establishing the Majorana nature 
of massive neutrinos at present are the (/3/3)o!y-decay experiments searching for the process 
{A, Z) — i> {A, Z + 2) + e~ + e~ (for reviews see, e.g. [131 HH [13 [I^)- The observation of 
(/3/3)o!/-decay and the measurement of the corresponding half-hfe with sufficient accuracy, 
would not only be a proof that the total lepton charge is not conserved, but might provide 
also unique information on the i) type of neutrino mass spectrum [T71 [18] (see also [191 [20]). 
ii) absolute scale of neutrino masses (see, e.g. [H]), and iii) Majorana CP- violating (CPV) 
phases [211 IISII221 [23] (see also the related discussions in, e.g. [SI [23 [2S] ) • 

Under the assumptions of 3-z/ mixing, of massive neutrinos uj being Majorana particles, 
and of (/?/5)oi/-decay generated only by the (V-A) charged current weak interaction via the 
exchange of the three Majorana neutrinos having masses rrij ^ few MeV, the {j3f3)oy- 
decay amplitude has the form (see, e.g. [131 [IB]): A{I3(3)qi, = <m> M, where M is the 
corresponding nuclear matrix element (NME) which does not depend on the neutrino mixing 
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parameters, and <m> is the (/9/3)oi.-decay effective Majorana mass, 

|<m>| = |mi |t/ei|2 + m2|t/e2pe^"^i +m3|[/e3|^e^"^^| . (2) 

Here \Uej\, j = 1, 2, 3, are the absolute values of the elements of the first raw of the PMNS 
mixing matrix, |f/ci| = C12C13, \Uc2\ = s^Cis, \Ue3\ = S13, Cij = cos%, Sij = sin%, 6112 = 9q, 
023 = Oa and 613 being the solar neutrino, atmospheric neutrino and CHOOZ mixing angles 
in the standard parametrisation of t/pMNS (see, e.g. [H]), and 0:21, 031 are the two Majorana 
CP- violation phases in t/pMNS [271 [28] . 

The experimental searches for (/3/3)o!/-decay have a long history [Hj. The best sensitivity 
was achieved in the Heidelberg- Moscow ''^Ge experiment [29]: |<m>| <(0.35 - 1.05) eV 
(90% C.L.), where a factor of 3 uncertainty in the relevant NME (see, e.g. [30]) is taken 
into account. The IGEX collaboration has obtained [31]: |<m>| < (0.33 — 1.35) eV 
(90% C.L.). A positive signal at > 3cr, corresponding to |<m>| = (0.1 — 0.9) eV, is 
claimed to be observed in [32j- Two experiments, NEM03 (with ^°°Mo and ^^Se) [33] and 
CUORICINO (with ^^OTe) [M], designed to reach a sensitivity to |<m>| ~ (0.2 - 0.3) 
eV, set the limits: |<m>| < (0.7 - 1.2) eV [33] and |<m>| < (0.2 - 0.9) eV P] (90% 
C.L.), where estimated uncertainties in the NME are accounted for. Most importantly, a 
large number of projects aim at a sensitivity to |<m>| ~(0.01-0.05) eV [35]: CUORE 
(^30Te), GERDA (^^Ge), SuperNEMO, EXO (^^exe), MAJORANA (^^Ge), MOON (^o^Mo), 
COBRA (ii^Cd), XMASS (^^^^Xe), CANDLES (^^Ca), etc. These experiments, in particular, 
will test the positive result claimed in [32] . 

The predicted value of | < m > | depends strongly on the type of z/— mass spectrum [TTl [18] , 
more precisely, on the type of hierarchy neutrino masses obey. Let us recall that the neutrino 
mass spectrum (in a standardly used convention) can be with normal ordering, mi < m2 < 
m3, or with inverted ordering, 7/13 < mi < m2. The first corresponds to Am\ = Am|-^ > 
0, |Am^| ~ (0.05)^ eV^ being the neutrino mass squared difference responsible for the 
(dominant) atmospheric neutrino oscillations; the second is realised if Am\ = < 0- 

Depending on the sgn(Am^) and the value of the lightest neutrino mass, i.e., the absolute 
neutrino mass scale, min(mj) = m^^j^, the neutrino mass spectrum can be 

i) Normal Hierarchical (NH): mi <^ m2 < m3, m2 = (AmQ)2, m3 = (Am^)2, AttIq = 
Am2i ~ 0.009 eV being the neutrino mass squared difference driving the solar oscillations; 

ii) Inverted Hierarchical (IH): m3 -C mi < m2, with mi 2 = |Am^|2, AtUq = Am^i; 

iii) Quasi-Degenerate (QD): mi = m2 = ms = mg, m| ^ |Am^|, mo ^ 0.10 eV. 

The existence of significant and robust lower bounds on | < m > | in the cases of IH and 
QD spectra [I?] (see also [IS]), given respectively by |<m>| ^ 0.01 eV and |<m>| ^ 0.03 
eV, which lie either partially (IH spectrum) or completely (QD spectrum) within the range of 
sensitivity of the next generation of (/3/?)o,y-decay experiments, is one of the most important 
features of the predictions of |<m>|. At the same time we have |<m>| ^ 5 x 10"'^ eV 
in the case of NH spectrum [23] • The fact that max{\<m>\ ) in the case of NH spectrum 
is considerably smaller than min{\<m>\ ) for the IH and QD spectrum opens the possi- 
bility of obtaining information about the type of z/-mass spectrum from a measurement of 
|<m>| 7^ [17]. More specifically, a positive result in the future generation of {I3(3)qu- 
decay experiments with |<m>| > 0.01 eV would imply that the NH spectrum is strongly 

^Up to small corrections we have in the cases of two spectra 17J: |<m>| ^ Am^cos26'o (IH) and 
[ < m > I ^ mo cos 26*© (QD) . The possibility of cos 26*0 = is ruled out at ^ 6(t by the existing data [36l [37] , 
which also imply that cos 26*0 ^ 0.26 at 2a [37]. We also have Am^ ^ 2.0 x 10"'^ eV^ at 3cr(see further). 
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disfavored (if not excluded). For Am\ > 0, such a result would mean that the neutrino mass 
spectrum is with normal ordering, but is not hierarchical. If Am\ < 0, the neutrino mass 
spectrum would be either IH or QD. 

If the future (/3/3)o,y-decay experiments show that |<m>| < 0.01 eV, both the IH and 
the QD spectrum will be ruled out for massive Majorana neutrinos. If in addition it is 
established in neutrino oscillation experiments that the neutrino mass spectrum is with 
inverted ordering, i.e. that Am\ < 0, one would be led to conclude that either the massive 
neutrinos Uj are Dirac fermions, or that I'j are Majorana particles but there are additional 
contributions to the (/5/3)ojy-decay amplitude which interfere distructively with that due to 
the exchange of light massive Majorana neutrinos. However, if Am\ is determined to be 
positive in neutrino oscillation experiments, the upper limit |<m>| < 0.01 eV would be 
perfectly compatible with massive Majorana neutrinos possessing NH mass spectrum, or 
mass spectrum with normal ordering but partial hierarchy, and the quest for |<m>| would 
still be open. 

If indeed in the next generation of (/3/3)oy-decay experiments it is found that |<m>| < 
0.01 eV, while the neutrino oscillation experiments show that Am\ > 0, the next frontier in 
the searches for (/?/3)o;y— decay would most probably correspond to values of |<m>| ~ 0.001 
eV. Taking |<m>| = 0.001 eV as a reference value, we investigate in the present article the 
conditions under which | < m > | in the case of neutrino mass spectrum with normal ordering 
would be guaranteed to satisfy | < m > | 0.001 eV. We consider the specific cases of i) normal 
hierarchical neutrino mass spectrum, ii) of relatively small value of the CHOOZ angle ^13 
as well as iii) the general case of spectrum with normal ordering, partial hierarchy and a 
value of close to the existing upper limit. We study the ranges of the lightest neutrino 
mass nil and/or of sin^ 613, for which |<m>| ^ 0.001 eV and discuss the phenomenological 
implications of such scenarios. 

In the present analysis we do not include the effect of the uncertainty related to the 
imprecise knowledge of the (/5/?)oi/— decay nuclear matrix elements (see, e.g. [30]). We hope 
that by the time it will become clear whether the searches for (/3/3)o^— decay will require a 
sensitivity to values of |<m>| < 0.01 eV, the problem of sufficiently precise calculation of 
the (/3/3)oi/— decay nuclear matrix elements will be resolved B 

The paper is organised as follows. In Section 2 we present predictions for |<m>| 
using the present 2a experimentally allowed ranges of values of the neutrino oscillation 
parameters and future prospective uncertainties in their values. In Section 3 we analise the 
conditions under which | < m > | in the case of u mass spectrum with normal ordering would 
be guaranteed to satisfy |<m>| ^ 0.001 eV. We consider the cases of i) normal hierarchical 
spectrum, ii) small ^13, and ii) spectrum with partial hierarchy. In Section 4 we give an 
estimate of | < m > | when the three u masses and the neutrino mixing originate from neutrino 
mass term of Majorana type for the (left-handed) flavour neutrinos and Yl^j'^j^lj = 0; but 
(/3/3)oy-decay is allowed. Section 5 contains the conclusions of the present analysis. 

•^Encouraging results, in what regards the problem of calculation of the NME, were reported in |30j . 
A possible test of the NME calculations is discussed in [38]. Let us note that nuclear matrix elements 
uncertainties do not affect the predictions for the effective Majorana mass parameter directly, but induce a 
spread on the values of the (/3/3)oi/-decay half-hfe times which correspond to the predicted values of |<77i>| . 
Conversely, if a measurement of the half-life time is performed or a stringent bound is obtained, they would 
affect the experimentally determined value of |<m>| and the constraints following from the latter. 
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2 Neutrino Oscillation Data and Predictions for | < m > | 

The existing neutrino oscillation data allow us to determine the parameters which drive the 
solar neutrino and the dominant atmospheric neutrino oscillations, Am^ = Am'^i, sin^^i2 = 
sin^6'0, and |Am|^| = |Am|;^| = lAm^gl, sin^26'23, with a relatively good precision, and to 
obtain rather stringent limits on the CHOOZ angle [39] ^13 (see, e.g. [36l|37j). The best fit 
values and the 2a allowed ranges of |Am^|, Am| and sin^ 9q read [TF] : 

(|Ami|)BF = 2.4 X 10"^ eV^ 2.1x10-3 eV ^ < |Ami| < 2.7 x lO^^ eV ^ (3) 
(Am| )bf = 7.6 X 10"^ eV ^ 7.3 x 10"^ eV ^ < Am| < 8.1 x 10"^ eV ^ (4) 
(sin^ ^0)bf = 0.32, 0.28 < sin^ Oq < 0.37 . (5) 

A combined 3-z/ oscillation analysis of the global neutrino oscillation data gives [3Zj 

sin^ < 0.033 (0.050) at 2a {3a) . (6) 

The existing data allow a determination of AttIq, siyl^ 9q and |Am^| at 3a with an error 
of approximately 8%, 22%, and 17%, respectively [37j. Future oscillation experiments will 
improve considerably the precision on these basic parameters: the indicated 3a errors could 
be reduced to 4%, 12% HI] and better than 5% SSI US] (see also the discussion in 
[m [23] and the references quoted therein), and even to ~ 1% for Am^[44j. "Near" future 
experiments with reactor i7e can improve the current sensitivity to the value of sin^ ^13 by a 
factor of (5-10) (see, e.g. [15]), while future long baseline experiments will aim at measuring 
values of sin^ 6*13 as small as 10~'^-10~3 (see, e.g. [HI US]). 

The type of neutrino mass hierarchy, i.e. sgn(Am^), can be determined by studying 
oscillations of neutrinos and antineutrinos, say, ^ z/g and ^ v^.) in which matter effects 
are sufficiently large. This can be done in long base-line z/-oscillation experiments (see, e.g. 
[IHl im 113]). If sin^ 26^13 ^ 0.05 and sin^ ^23 ^ 0.50, information on sgn(Am3i) might be 
obtained in atmospheric neutrino experiments by investigating the effects of the subdominant 
transitions z/^(e) v^^^^ and v^(^e) ^'e(/i) of atmospheric neutrinos which traverse the Earth 
[47] . For v^{e) {or i^^(e)) crossing the Earth core, new type of resonance-like enhancement of 
the indicated transitions takes place due to the (Earth) mantle-core constructive interference 
effect (neutrino oscillation length resonance (NOLR)) [48J 0. For Am^^ > 0, the neutrino 
transitions z/^(e) z/e(^) are enhanced, while for Amg^^ < the enhancement of antineutrino 
transitions i/^(e) — >■ z/e(^) takes place, which might allow to determine sgn(Am|^). If sin^ ^13 
is sufficiently large, the sign of Am\ can also be determined by studying the oscillations 
of reactor u^. on distances of ~ (20 — 40) km ^T\. An experiment with reactor z/g, which, 
in particular, might have the the capabilities to measure sgn(Am^), was proposed recently 
in [52]. According to [52], this experiment can provide a determination of |Am^| with an 
uncertainty of (3 — 4)% at 3a. 

As is well-known, neutrino oscillations are not sensitive to the absolute scale of neutrino 
masses. Information on the absolute neutrino mass scale can be derived in /3-decay 
experiments [53t [5^ [55] and from cosmological and astrophysical data. The most stringent 

*As a consequence of this effect the corresponding ^^^(e) {or i^p(e)) transition probabihties can be maximal 
[49] (for the precise conditions of the mantle-core (NOLR) enhancement see [48l |49]). Let us note that 
the Earth mantle-core (NOLR) enhancement of neutrino transitions differs [48] from the MSW one. It also 
differs [48l from the parametric resonance mechanisms of enhancement discussed in the articles [50] • 
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upper bounds on the z/g mass were obtained in the Troitzk [M] and Mainz [55] experiments: 



rrir 



< 2.3 eV at 95% C.L. 



(7) 



We have rrip^ = mi^2,3 in the case of the QD u-mass spectrum. The KATRIN experiment [55] 
is planned to reach a sensitivity of rrip^ ~ 0.20 eV, i.e. it will probe the region of the QD 
spectrum. Information on the type of neutrino mass spectrum can also be obtained in (3- 
decay experiments having a sensitivity to neutrino masses [56] ~ v^lAmil = 5 X lO'^ eV 
(i.e. by a factor of ~ 4 better sensitivity than KATRIN [55]). 

The CMB data of the WMAP experiment [57], combined with data from large scale 
structure surveys (2dFGRS, SDSS), lead to the following upper limit on the sum of neutrino 
masses (see, e.g. [58]): 



S < (0.4-1.7) eV at 95% C.L. 



Data on weak lensing of galaxies, combined with data from the WMAP and PLANCK 
experiments, may allow S to be determined with an uncertainty of ~ 0.04 eV 
It proves convenient to express 



the three neutrino masses in terms of Am'i, and 







Ami, measured in neutrino oscillation experiments, and the absolute neutrino mass scale 
determined by min(mj) = mj^jp^lf]. In both cases of z/-mass spectrum with normal and inverted 



ordering one has (in the convention we use) 
normal ordering, m^jj^ = mi, 



AuIq 



Am|i > 



Ami 



Amg^ > and ms 



m2 

ml + Ami 



ml + Am| 



spectrum is with inverted ordering, m^jj^ = m^, Ami 



Am32 < and mi 



For 

5, while if the 
(mo + lAmA I — 



Am; 



2 U 



For the elements of the PMNS matrix 



\U -P 



1,2,3, as we have already 



indicated, the following relations hold: |t/cip = cos^ 6^0(1 — sin^ 6* 
sin^6'i3), and |?7e3p = sin^6'i3. Thus, given |Ami|, Am^, 9q and 9 



13) 



|<m>| 



sin2^0(l - 
713, depends on 

the lightest neutrino mass (absolute neutrino mass scale), m^jj^, the two Majorana phases 
Q!2i and ^31, present in the PMNS matrix and on the type of neutrino mass spectrum (see, 
e.g. ^18j). For neutrino mass spectrum with normal ordering we have 



<m> 



m^j^cos 



;i - sin^ ^13) + Jml^+ Ami 



sin^ 9Qi 



'9 



13 J 



+ A/m,2^,+ Amisin2^i3e*"^i 



(9) 



For spectrum with inverted ordering a different expression is valid [211 flS 



\<m>\ 



m^jj^+ I Ami I ~ Am^ cos^ 6'0(1 — sin^ 9 



13 j 



+ \h^mm+ I Amil sin^ ^0(1 - sin^ ^i3)e'"^^ + m^^^sin^ ^136 



(10) 



m^j^+ I Ami I |cos^ 9q + sin^ ^0 e^"^^ 1(1- sin^ ^13) , m^^^ = m3 . (11) 



^For a detailed discussion of the relevant formalism see, e.g. [T51 116j . 
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In Eq. (ITT]) we have neglected Am| with respect to (m^^^ + |Am^|) and the term 
m^jj^sin^ 6'i3. According to the existing data, we have AttIq /(^«min + I^^^aI) ^ 0.032, 
and m^jj^sin^ ^^13 ^ ("^min+ I^^^aI)^ cos 26^0. Actually, the term m^jj^sin^ 6'i3 can always be 
neglected provided sin^ ^13 -C cos 29q. The expression for | < m > | in the case of IH spectrum 
follows from Eq. ffTTj) if m^jj^-C |Am^| and m^jj^ is neglected with respect to |Am^|. For 
the QD spectrum we get: 

|<m>| = mo Kcos^^Q + sin^^oe*"^!) (l-sin2ei3) + sin2ei3e^°^i| , (12) 

^ mo Icos^ Oq + sin^ Oq e*"^^ | (1 - sin^ ^13) , (13) 

where mo = "^^1^ "^i — "^2 — f^s- Evidently, as long as sin^ ^^13 ^ cos2^^0, the terms 
oc sin^ ^13 in | < m > | play an insignificant role in the cases of neutrino mass spectrum with 
inverted ordering (i.e. Am^ < 0), or of QD type (for any sgn(Am^)). In what concerns 
the spectrum with normal ordering, the term \/ ^^1^™+ Am^ sin^ 6*13 can be crucial for de- 
termining the magnitude of |<m>| if massive neutrinos are not QD, i.e. if m^jj^^ Am^, 
and sin^ 6'i3 is sufficiently large (see further). 

If CP-invariance holds, we have [61] 021 = kn and 031 = k'n, k,k' = 0,1,2,.... In the 
case of CP-invariance the phase factors 

7^21 = e'°2i = ±1 , r/31 = e^"3i = ±1 , r];^ = e'"^^ = ±1 , (14) 

as is well-known, have a simple physical interpretation [6II |T3]: r]ik is the relative CP-parity 
of Majorana neutrinos z/j and z/^. Obviously, |<m>| depends strongly on the Majorana 
CPV phase(s): the CP-conserving values of 021 = 0, ±7r determine, for instance, the range 
of possible values of | < m > | in the cases of IH and QD spectrum. 

We recall that the neutrino oscillation experiments are insensitive to the two Majorana 
CP- violation phases in the PMNS matrix f27\ [52] - the latter do not enter into the expressions 
for the probabilities of flavour neutrino oscillations. It is interesting to note, however, that 
in addition of playing an important role in the predictions for |<m>| and, correspondingly, 
of the (/5/5)oi/— decay half-life, the Majorana phase(s) in f/pMNS can provide the CP-violation 
necessary for the generation of the baryon asymmetry of the Universe [631IM] (see also [65j). 
The Majorana phases 021 and 032 can also affect signiflcantly the predictions for the rates 
of (LEV) decays yU^e-l-7, r^/i-l-7, etc. in a large class of supersymmetric theories with 
see-saw mechanism of u-mass generation [66] . 

Eirst, we will update the predictions for |<m>| as a function of m^jj^, using as input the 
2a ranges of values of Am^, Am^ , sin^ ^0 and sin^6'i3, obtained from the latest available 
set of neutrino oscillation data (see Eqs. ([3]), (jl]) and (^). Since a2i and 031 cannot be 
determined in independent experiments, we treat them as free parameters taking values 
< 0:21,31 < 271. The results of this analysis are shown in Eig. 1. 

We report in Table 1 the maximal and minimal values of |<m>| for the normal hi- 
erarchical (NH) spectrum, mi -C m2 < m3, for the inverted hierarchical (IH) spectrum, 
m3 <^ mi ~ m2, and for the quasi-degenerate spectrum (QD), mi ~ m2 ~ m3 > 0.2 eV. 

In Eig. [2] we show the predicted ranges of |<m>| using the present best flt values of the 
neutrino oscillation parameters and their prospective errors as discussed above. We assumed 
a la experimental error of 2%, 2% and 4% on AmQ, Am^, sin^ 9q, respectively. Eor sin^ 6*13, 
we take sin^ ^13 = 0.01 and we consider the la uncertainty in the absolute value of 0.006. In 
Table 2 we give the maximal and minimal values of |<m>| for the three spectra, NH, IH 
and QD. 
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1 ^ ^ 1 NH 
l<"^>l min 


1 ^ ^ 1 NH 

\<m>\ 


\<m>\ 


l<"^>l max 


l<m>| 

1^ "'^1 min 


0.7 


4.8 


11.3 


51.5 


44.2 



Table 1: The maximal values of |<m>| (in units of meV) for the NH and IH spectra, and 
the minimal values of |<m>| (in units of meV) for the NH, IH and QD spectra, obtained 
using the 2a allowed values of the neutrino oscillation parameters. The results for the NH 
and IH spectra are for m^^-^ = 10~^ eV, while those for the QD spectrum correspond to 



1 ^ ^ 1 NH 


1 ^ ^ 1 NH 


1 ^ ^ 1 IH 
l<"^>l min 


1 ^ ^ 1 IH 
l<"^>l max 


l<m>| 

1^ '"^1 min 


2.1 [1.5] (1.0) 


3.5 [3.9] (4.4) 


15.1 [15.0] (14.8) 


50.1 [50.0] (49.1) 


63.4 [60.7] (58.0) 



Table 2: The maximal values of |<m>| (in units of meV) for the NH and IH spectra, and 
the minimal values of |<m>| (in units of meV) for the NH, IH and QD spectra, for the best 
fit values of the oscillation parameters and using the prospective errors discussed in the text. 
We take sin^ ^13 = 0.0 [0.01] (0.02). The results for the NH and IH spectra are obtained for 
'^MiN" while those for the QD spectrum correspond to rn^^^ = 0.2 eV. 



3 The |<m>| ^ 10 ^ eV Frontier in (/?/?)ozy— Decay 

In the present Section we will analise the conditions under which |<m>| ^ 10"'^ eV in the 
case of neutrino mass spectrum with normal ordering. Before discussing the general case of 
arbitray mi and sin^ ^13 satisfying the presently exiting experimental limits, we will consider 
two specific but physically interesting cases: i) negligibly small mi (NH spectrum), and 
ii) relatively small sin^^i3, such that the term a/ mf + Am\ sin^ 6'i3 in Eq. ([9]) is strongly 
suppressed, i^mf + Am\ sin^ 6'i3 ^ 10~^ eV. 

3.1 Normal Hierarchical Spectrum 

In the case of the normal hierarchical spectrum we have mi -C m2,3 and therefore only the 
two heavier neutrinos, z/2 and 1/3, contribute to the effective Majorana mass parameter. In 
this case m2 = ^^/Am^, = -^/AmJ, and the sum of neutrino masses reads: 

mi + ma + m3 = 0.058 eV . (15) 

The effective Majorana mass is given by: 



|<m> 

where 032 = 031— a2i is the difference of the two Majorana CP- violating phases in f/pMNS- We 
will refer to the first term in the r.h.s. of Eq. (fT6|) as the "solar term" due to its dependence 
on AttIq , while to the second as the "atmospheric" one. The two terms in the expression for 
|<m>| add constructively if < 032 < 7r/2, while for tc/2 < 032 < n partial or complete 
cancellation between the "solar" and "atmospheric" terms can take place. The cancellation 



Am| 



sm 



^0(1 



sm' 



'^13) + 



sm 



iaz2 



(16) 
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is most effective in the case of CP-invar iance and 032 = vr. The degree of cancellation is 
controlled by sin^6'i3. For sufficiently small values of 61^, sin^ 6'i3 ^ 0.01, the solar term 
dominates and |<m>| is predicted to be in the few meV range, |<m>| ~ (2 — 3) x 10^'^ 
eV. If sin^6'i3 is close to the present 3a bound [37], sin^ 6'i3 < 0.05, the solar and the 
atmospheric terms in Eq. ( |T6l) are of the same order and a substantial cancellation can take 
place. We will analise this possibility first qualitatively. 

Consider the "extreme" case of 032 = vr and | < m > | = [§. This requires [HI [191 126] 

|<m>|=0: sin^ ^13 = V^!!^ gin^ Oq , (17) 

where we have neglected AitLq with respect to Am^. Taking the best fit values of Ami, , 
sin^ 9q and Am^, determined from the analysis of the currently existing neutrino oscillation 
data, we get sin^ 613 = 0.057, which is ruled out by the data. Using the 2a and 3a ranges 
of allowed values of the same three parameters, we find respectively sin^ ^13 = 0.046, which 
is close to the current 3a upper limit on sin^6'i3, and sin^ 6*13 = 0.041. Thus, in order 
for |<m>| to be strongly suppressed, |<m>| <^ 10~^ eV, sin^ 6*13 should have a value 
close to the existing 3a upper limit. If we use the current 2a {3a) upper limit on sin^^i3, 
sin^ 6'i3 < 0.033 (0.050), and the present best fit values of AuIq , siia^ 9q and Am\, we 
find for 0:32 = tt that |<m>| ^ 1.1 (0.2) x 10"^ eV. If < 0:32 < 57r/6, we obtain 
|<m>| ^ 1.5 (1.3) X 10~^ eV. It follows from this simple analysis that if, in the future high 
precision measurements of AttIq, sin^ ^0 and Am\, the currently determined best fit values 
of these parameters will not change and sin^ 6*13 is found to have a value sin^ 6'i3 ^ 0.01 (0.03), 
the effective Majorana mass will satisfy |<m>| ^ 2.2 (1.2) x 10^'^ eV for any 032. For, 
e.g. < 0^32 < 57r/6, we have |<m>| ^ 1.3 x 10^^ eV for any sin^ ^13 allowed at 3a by the 
existing data. Values of 032 7^ in the indicated range are required for the generation of 
the baryon asymmetry of the Universe in the "flavoured" leptogenesis scenario, in which the 
requisite CP- violation is provided exclusively by the Majorana phase (difference) 032 [63]. 

We will perform next a similar analysis of the conditions under which |<m>| ^ 10^'^ 
eV, taking into account the current and prospective uncertainties in the measured values 
of the relevant neutrino oscillation parameters. The minimal predicted value of |<m>| , 
|<m>|^^^ , is obtained in the case of CP-conservation and opposite CP-parities of the two 
relevant neutrinos and can be evaluated as 

|<m>| = |<m>| —na(\<m>\ ), (18) 

where | < m > | _ is the predicted value of | < m > | obtained using the best fit values of the 
oscillation parameters, a{\<m>\ _) is the error on |<m>| and n = 1,2,3.... 

Using the propagation of errors and assuming that the errors on the oscillation parameters 



^We postpone the discussion of the decay in the case of |<m>| = to Section 4. 
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of interest are small and independent, we obtain the 1-a error on |<m>| for any 032: 

a{\<m>\) ~— — ^ 1 [ sin'' 6113 Am\ ( sin^ 6 13 J Am| + J Am| sin^ 6q cos 032] (5^(Am|) 

+ Am| sin^ Oq (^^/Am^sm^ Oq + A/Am^sin^ 6*13 cosa32)^ (^452(31^2 ^^^^ _^ 52(^Am| ) j 

+4 f sin^ 6^13 Am^- Am| sin'' ^0 + A/Am| Am| sin^ cos 032 ) (T^(sin^ 6'i3)j 

(19) 

Here 5(sin^ 6'0), 5(AmQ) and 6{Am\) are the relative errors on the oscillation parameters 
Amg, sin^ 6*0 and Am\, a^siia^ 613) is the absolute error on sin^^i3, and we have used the 
fact that sin^ 6*13 -C 1. We have assumed (see Section 2 and Fig. 2) and will use in our further 
analysis (see Section 2) the following values of the errors: 5(sin^ 9q) = 4%, 6{AmQ) = 2% 
and S{Am\) = 2%. For the chosen values 5(sin^ ^0), 5(AmQ) and 6{Am\), the error on 
AitLq gives a sub dominant contribution in comparison with that on the solar mixing angle 
and we neglect it in the following discussion. 

If CP-invariance holds we have 0:32 = 0, vr and Eq. (fT9l) simplifies to: 



sin 6^'iAm?-, 



a(\<m>\ ^) ~ \ Ami sinH q5\siti^ 6 q) + ^ 5\Aml) + Amia^(sin^ ^13) 



4 

(20) 

where we have neglected -^/Am^sin^^© with respect to a/ Am\. In Eq. (!20l) \<m> 



± 



refers to 7732 = ±1. The contribution of the error on Am^ in a{\<m>\ ^ is suppressed 



by the factor sin^ ^^13 and can also be neglected, while the errors on sin^ 613 and on sin^ 6q 



±> 

and on sin^ 

can give sizable contributions to a{\<m>\ ^ and both should be taken into account. For 
the current best fit values of the oscillation parameters, a{\<m>\ ^) is given to a good 
approximation by cr(|<m>| ^ ) ^ v^A^v^(0.057(52(sin" ^0)) + cr2(sin ^13). It is clear from 
this expression that for an error on sin^ ^0 of 4-8%, the two terms in cr(|<m>| ^) are of the 
same order if cr(sin^ ^^13) = 0.004, while for cr(sin^ ^^13) 0.006 the error on sin^ ^13 typically 
gives the dominant contribution in a{\<m>\ _|_). 

For neutrinos of equal CP-parities, i.e. ^32 = 0, the mean value of |<m>| is predicted 
to be in the few meV range and the expected relative error ij{\<m>\ _|_) varies between 7% 
and 15%, depending on the specific values of errors and best fit values of the parameters. If 
the neutrinos z/2 and z/3 have opposite CP-parities, i.e. 0:32 = vr, the mean value of |<m>| is 
smaller as partial cancellation between their contributions to | < m > | can take place. In this 
case the error on |<m>| can become as large as 30%-40%. 

If CP-symmetry is broken, the full expression for a{\<m>\ ), Eq. f[T^ . should be used. 

It can be shown, however, that a{\<m>\ ) < max^cr(|<m>| ^),cr(|<m>| _)^. 

Using Eq. fl20|) in the case of 7732 = — 1, we can study analytically the condition on sin^ 613 
which guarantees that the predicted value of |<m>| is larger than 1 meV. Neglecting the 
dependence on sin^ ^13 in cr ( | < m > | ) , we find an approximate solution for sin^ ^13 : 



A/Am| sin^ 6*0-1 meV - n^^| Ami ^^^"^ OQd'^ism^ 6q) + Am|cr2(sin^ 613) 
sin^ 6*13 < 



V^Am^ 



A 

(21) 



10 



In Fiff. 3 we show the values of sin^^T? versus Ami for which |<m>| =1 meV 

to -LJ A I I MIN 

is satisfied for n=l, 2, 3 (dash-dotted, dashed, dash-double-dotted lines). We use the best 
fit value of sin^ 6'© and two values of the error on sin^^is. If sin^ ^13 is larger than the 
shown values, a strong cancellation between the two contributions to | < m > | can take place 
and |<m>|^^^ < 1 meV. This would imply that, depending on the value of 032, there 
are predicted values of | < m > | both smaller and larger than the future reference sensitivity 
used in this analysis. The possibility for a future experiment to find a positive signal of 
(/3/5)ov-decay would depend on the unknown value of 032. 

The limiting value of sin^ 613 is in the 0.01-0.03 range. The precise value depends critically 
on the error on sin 61^: for cr(sin ^13) ~ 0.004 (0.008), we have sin^ ^13 < 0.02 (0.01). The 
limit on sin^ 613 depends also on Am\, as can be easily understood from Eq. (12T1) : the larger 
Am\, the smaller the bound on sin^6'i3. The value of sin^ 9q controls the magnitude of 
the first term in |<m>| _ and therefore plays an important role in Eq. fl2Tl) . We show the 
dependence on sin^ 6q in Figs. 4 and 5. The smaller the value of sin^ 6q, the smaller sin^ ^13 
for which one can have |<m>| < 1 meV. If sin^ Op, = 0.26 and cr(sin^^i3) ~ 0.004, we 

I I MIN ^ \ -LO/ 1 

have |<m>| > 0.001 eV for values of sin^6'i3 < 0.01. If, however, a{sm^ 9is) — 0.008, one 
can have |<m>|^^^ < 1 meV even if the (mean) value of sin^ 6'i3 = 0. On the contrary, 
for sin^ 9q = 0.40, a large part of the relevant parameter space is already excluded by the 
present data [37] and we get |<m>| > 0.001 eV for sin^ 6'i3 < 0.03 (0.02) in the case of 
a{sm^9i3) = 0.004 (0.008). 

The preceding rather detailed analysis shows that |<m>| > 0.001 eV typically for 
sin2 6'i3 ^ (0.01-0.02). Values of sin2 6'i3 (0.01-0.02) are within the sensitivity of the two 
reactor experiments Double-CHOOZ [67j and Daya Bay ^j, which are under preparation, 
and of the currently operating and future long baseline neutrino oscillation experiments, 
MINOS [8], OPERA [29], T2K and NOi/A [16]. The results of these experiments will be 
crucial for establishing whether the effective Majorana mass |<m>| in the case of NH 
neutrino mass spectrum is limited from below and for determining its lower limit. 

In the case of a NH spectrum, for sin^^is = 0, only one contribution in |<m>| is 
relevant, the other two being suppressed by the negligible values of mi and sin^ 6*13. In this 
case there is no dependence of |<m>| on 0:32. If sin^ ^13 has a value close to the existing 
upper limit, a sufficiently accurate measurement of |<m>| could allow to distinguish the 
two possible CP-parity patterns or establish CP-violation. Here, we study what would be the 
requirements in order to have sensitivity to CP-violation. We perform a simplified analysis 
in which we retain for both CP-parity patters only the dominant term in the theoretical 
error on | < m > | : 



a(|<m>| ) ~ yAm|(T(sin^^i3) . (22) 

The existence of a "just-CP-violating" region [18], signaling the possibility to search for CP- 
violation, requires the allowed regions for the CP-conserving cases 7732 = 1 and r/32 = — 1 not 
to overlap. This condition is satisfied provided 

sin^^i3 > na(sin2^i3), (23) 

where n is the number of ai^sm"^ 613) considered. For example, for a^sin^ 613) = 0.004,0.008 
and n = 2, we have sin^ 6*13 > 0.008, 0.016. In this case, in principle, it would be possible to 
distinguish the two CP-parities patterns or find CP-violation due to a Majorana CP- violating 
phase. CP- violation would be established if the experimentally allowed value of |<m>| is 
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within the "just-CP- violating" region, once the experimental error on |<m>|, A, and the 
nuclear matrix elements uncertainties are taken into account. Even if Eq. fl23l) is satisfied, 
this is a formidably challenging task. In the most optimistic case of sin^ Oi^ having a value 
close to the present 3cr bound, sin^ ^13 ^ 0.05, for a nuclear matrix element uncertainty 
( = 1.5 on |<m>|, an error not larger than A = 0.5 meV would be required. The width of 
the "just-CP- violating" region decreases rapidly with ^13 and for smaller values of sin^ 613 
the error required on | < m > | would be even smaller. 



3.2 The Case of Small sin^ ^13 

Consider next the possibility of sin^ ^13 having a rather small value, such that 
y^m|~TAm|" sin^ 6^13 < 2 x 10~^ eV -C 10"^ eV. For ml <^ Am\ this conditions is fulfilled 
if sin^ ^13 ;^ 4 X 10"^, while if, e.g. mi = 0.05 eV, it is satisfied provided sin^ 6'i3 ^ 3 x 10"^. 
These values of sin^ 613 can be tested, e.g. in future long baseline neutrino experiments with 
superbeams, beta beams and at neutrino factories |1T[ HS] . 

We set sin^ ^13 = for simplicity in the following discussion. The expression for |<m>| 
simplifies to: 

|<?7i>| = mi cos^6'm + A/mf+Am^ sin^6'me*""i . (24) 



For NH neutrino mass spectrum, i.e. for mi <^ ^jAm^ sin^ 6q/ cos^ 6q = 4 x 10~^ eV, we 

4 



always have |<m>| = 3x10^ eV. If, however, ml ^ Am^ , the neutrino mass spectrum 



will not be hierarchical. There are two possibilities, 
i) For ml > Am| = 7.6 x 10"^ eV^, we get 



|<m>| = mi y 1 - sin^ 26'0 sin^ ^ ^ mi cos 26^0 . (25) 

Taking mi 2 x 10"^ eV and the 2a {3a) lower limit o n cos26'q, cos 26)0 > 0.26 (0.20), we 
find |<m>| ^ 5.2 (4.0) x 10~^ eV. In this case m2 = ^ml + Am| ^ 2.2 x lO'^ eV, and 
the sum of neutrino masses satisfies: 

mi + m2 + m3 9.5 X 10~^ eV . (26) 

ii) If, however, ml ~ Amg and 021 ~ vr, a cancellation between the two terms in Eq. 
(!24|) is possible and |<m>| can be strongly suppressed, |<m>| <^ 10~^ eV. Consider the 
extreme case of |<m>| =0 (for a more detailed discussion of the (/3/5)o;/— decay in the case 
of |<m>| =0 see Section 4). For 021 = tt, it is reahsed if [13 

|<m>| =0: mi = m2 tan^ 6*0 . (27) 



Using the relation m2 = {ml + Am^ ) 2 , we find that | < m > | =0 can hold in the case being 



1 

''0 ' 

studied if mi = Amg sin^6'0/ cos 2^0 = 2.2 x 10~^ eV^, where we have used the best fit 
values of Am^ and sin^6'0. This implies that mi = 4.6 x 10~^ eV, m2 = 10~^ eV, and, 
corresp ondingly, 

mi + m2 + m3 = 6.4 X 10"^ eV. (28) 
It is not difficult to convince oneself, however, that if 021 = tt, one obtains |<m>| ^ n for 



mi ^ 



cos 29(r 







cos 9q + siir 9q + ji-'^ Am| cos 2610 



(29) 
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where the reference value /i = 10~^ eV in the case of interest. Using the best fit values of 
and sin^ 6q we get mi ^ 6.6 x 10~^ eV. For the sum of neutrino masses we obtain 
mi + m2 + m3 ^ 6.7 x 10~^ eV. 

This qualitative analysis shows that if sin^ ^13 ^ 3 x 10^^ and the sum of neutrino masses 
satisfies mi + m2 + ms 7 x 10^^ eV, we will have |<m>| > 10^^ eV for any 021 • 



3.3 Spectrum with Partial Hierarchy 



As is well-known |T9], in the case of neutrino mass spectrum with partial hierarchy we can 
have |<m>| ^1 meV and even |<m>| =0. However, this requires that the lightest 
neutrino mass mi has a value in the rather narrow interval, mi ~ (few x 10"'^ — 10~^) eV. As 
a consequence, the sum of neutrino masses should also lie within a specific interval. Here we 
analyze the values of mi and sin^ ^13 for which the indicated strong cancellation in | < m > | 
would not take place and we would have | < m > | > 1 meV. 

For the neutrino mass spectrum under discussion, all the three contributions to | < m > | in 
Eq. (ini) are relevant. We consider the effect of cancellations between the three terms in the 
case of CP-invariance, in which there are four different neutrino CP-parity patterns. We will 

denote them as + + + (H ) if 0:21,31 = (vr), and + H — (H h) when 0:21 = (vr) 

while 031 = TT (0). The prediction in the case of CP- violation will lie within the ones 
obtained for CP-conservation. Obviously, if both ^ 0:21 ^ vr/2 and ^ 0:31 ^ 7r/2, there 
will be no mutual compensation between the three terms in Eq. ([9]) and we would have 
|<m>| 3 X 10-3 eV. 

For each CP-parity pattern, we analyse what are the values of mi and sin^ 6'i3 which would 
guarantee | < m > | > 1 meV or, conversely, which would be implied by a negative result for 
a search of neutrinoless double beta decay with a sensitivity of 1 meV, in the hypothesis 
of Majorana neutrinos. The effective Majorana mass parameter would be predicted to be 
smaller than 1 meV, if a sufficient cancellation between the three terms in the r.h.s. of 
Eq. ^ takes place. 

Here we use = 1 meV as a reference value for |<m>| , but similar results can be 
obtained for other values of jj, in the few meV range 0. The central value of mi can be found 
by solving Eq. ([9]) with | < m > | =1 meV, while the error on mi is obtained by propagating 
the errors on the oscillation parameters: 

f d\<m>\ , , a{\<m>\ ) , . 

^ mi = ^ a <m> ^ ^ . 30 

V 5mi ; ^^^^^^^ mism^e 



& 



The degree of cancellation between the three terms in |<m>| depends on the neutrino 
CP-parity pattern. The results for mi for the different CP-parity patterns are presented in 
Fig. 6 for three values of ^0, sin^^© = 0.26,0.32,0.40, using the prospective relative errors 
of 2%, 2% and 4% for Am\, AtUq and sin^ 9q and the absolute error of 0.006 on sin^ 6*13. 



We can understand the results in Fig. 6 by performing a simplified analysis neglecting 
cr(mi). We study each CP-parity pattern separately. In the following, we will use the present 
best fit values of Am^ , Am^ and sin^6'Q, unless otherwise indicated. 

^Let us note that a similar analysis for /i = was performed in Ref. [19] . 
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For the CP-parity pattern (+ + +), no cancellation takes place and we will have 
|<m>| ^ 2.5 meV for any allowed value of 613 and 6q. A negative search for 
neutrinoless double beta decay with a sensitivity of few meV, such that |<m>| q < 
|<m>| _,_ — ncr{\<m>\ ), where |<m>| q is the experimentally determined value of 
|<m>| and |<m>| ^ corresponds to NH spectrum and 7732 
would strongly disfavour (if not rule out) this possibility. 



+1 (see Section 3.1), 



In the (+ - 
term, a/A 



m. 



gnificant cancellation can take place only if the atmospheric 
sin^^is, is of the same order as the sum of the first two terms in the 



r.h.s. of Eq. ([9]). We can have |<m>| > jj =1 meV for given mi provided sin^ 6' 
satisfies 

mi cos^ 6*0 + ^/ml + AtUq sin^ 9q 
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sin^ 9i3 < 



(31) 



A/m^ + Am\ 

The above inequality is always fulfilled for mi ^ ^/Amg . For mi <^ A/Am^ , this 
condition becomes sin^ 613 < (sin^ ^^13)0 with 



(sin^ e 



no 



^ ■ sm^ 9p 







(32) 



It is satisfied for the best fit values of Am^, sin^ and Am^ , while if one uses the 3a 
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allowed ranges of these parameters, the inequahty implies sin^ 6*13 ^ 0.026, 0.036, 0.051 
for sin^6'0 = 0.26,0.32,0.40. These values of sin^6'i3 are close to the present 3a upper 
bound. In summary, for values of sin^6'i3 < (sin^6'i3)o, |<m>| is guaranteed to be 
larger than 1 meV for any mi. 

2fl ^ ^r.;^'iQ ^ have |<m>| > u =1 meV if mi satisfies 



For given sin 6*13 > (sin 6*13)0, 
mi > (m^)_, where 



m 



1 )± 



cos 26^0 



Am^ sin^ ^13 + j cos^6'0 



± 



sin^6'0 



{^^J Am\ sin^6'i3 + /i j + Am| cos26'0 



(33) 



In deriving Eq 

2/3 ^ \_ The lower bound 
10-3 



cos C713 

rather small: for /i 



331) we have neglected m{ with respect to Amj^ and have taken 
mf)_ of mi in Eq. fl33l) increases with sin^6^i3, but is 
^ — 0.05 and best fit values of the other relevant 



1 

eV, sin^ 6'i3 



oscillation parameters we get 



m 



1 )- 



0.9 X 10-3 eV. 



For the CP-parity pattern (H h), a partial cancellation can take place between 

the first and the second terms in Eq. ([9]); the cancellation would be significant only 
if mi ~ few meV. The second term in Eq. (|9]) would dominate and we would have 
|<m>| > yU = 10-3 eV only if mi and sin^ 6*13 are sufficiently small, more precisely. 



if < mi < 



mi 



and sin 6^13 < (sin 6*13)0, where (mf )_ and (sin 6*13)0 are given 



in Eqs. ( |33!) and (132!) . respectively. 
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The sum of the first and third terms in Eq. Q will dominate and will lead to 
|<m>| > fi = 10~^ eV for mi > (mf )_|_, where 



— \r- \( - \A"^A sin^ 6113 ) cos^6'0 



± sin^6'0 - y^Am^ sin^ 6113 j + Am| cos26'0 , (34) 



provided 



^13 < sin^ + . (35) 



Ami v/A< 



Given the experimental ?>a upper bound sin^ 6*13 < 0.05, the second inequality is always 
satisfied for /i = lO'^ eV. For sin^ ^13 = (0.02) we get from Eq. rrii^ 6.6 (4.7) x 
10-3 eV. 

Finally, consider the case (H ). As the second and third term in the r.h.s. of 

Eq. ^ are summed constructively, a strong cancellation in |<m>| can happen only 
for sufficiently large values of mi. We get |<m>| > /i = 10"^ eV for < mi < 
(— (mf )_) and for mi > (m^)_|_, where (m^)+ is given in Eq. fl33|) . The maximal 
value of mi determined by Eq. flMl) can be rather large. More specifically, we have 
-(mf )_ = 2.8 (5.0) [7.6] x 10"^ eV for sin^ ^13 = (0.025) [0.05]. For the minimal value 
of mi determined by the inequality mi > (m^) + , we get for sin^^i3 = (0.025) [0.05]: 
(m^)+ = 6.6 (9.3) [12.1] x 10^'^ eV. In the latter case the sum of the neutrino masses 
is limited from below by (mi + m2 + m3) ^ 6.8 (7.2) [7.9] x 10~^ eV. Both (— (m 
and (m^)+ increase with ^13 and sin^6'0. 



It follows from the preceding discussion that if a future highly sensitive (/3/3)oi/-decay 
experiment does not find a positive signal down to |<m>| ~ 1 meV, Majorana neutrinos 
would still be allowed, but the spectrum would be constrained to be with normal ordering 
and mi would be bound to be smaller than ~ 10~^ eV. The CP-parity pattern (+++) will 
be strongly disfavored (if not ruled out) as well. If in addition it is found that sin^ ^13 ^ 0.01, 
i) the CP-parity pattern (+ H — ) will also be disfavored, and ii) mi would be constrained 
to lie in the interval mi ~ (10^^ — 10~^) eV. No other future neutrino experiment will have 
the capability of constraining the lightest neutrino mass (and the absolute neutrino mass 
scale) in the meV range. Obviously, the above limits would hold only if massive neutrinos are 
Majorana particles. If the lightest neutrino has a mass in the interval mi ~ (10"'^ — 10~^) eV, 
this can have important effects on the generation of the baryon asymmetry of the Universe 
in the "flavoured" leptogenesis scenario of matter-antimatter asymmetry generation 



4 (/?/3)oz/-Decay in the Case of |<m>| =0 

In the present Section we shall discuss briefly the possible implications of having | < m > | = 
for the process of (/3/3)oi/— decay. If |<m>| = as a consequence of conservation of 
certain lepton charge, which could be, e.g. L^, L, or V = L^, — — L^., the (/?/?)o;y— decay 
will be strictly forbidden. However, in the case of neutrino mass spectrum with normal 
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ordering, one can have |<m>| = 0, as we have seen, as a consequence of an "accidental" 
relation involving the neutrino masses, the solar neutrino and CHOOZ mixing angles and 
the Majorana phase(s) in f/pMNS- For the spectrum of the normal hierarchical type, the 
relation of interest is given by Eq. i^, while if sin^^ig is negligibly small it is shown in 
Eq. ( 1271) . None of the two relations can be directly associated with a symmetry which 
forbids (/?/?)o;/— decay. Thus, if |<m>| = is a consequence of Eq. ( fT7|) or Eq. ( 1271) . 
(/?/5)oiy— decay will still be allowed. In what follows we will estimate the non-zero contribution 
to the (/5/3)ov— decay amplitude A{PP)qi, due to the exchange of the light massive Majorana 
neutrinos uj in the case when | < m > | =0 and there is no symmetry forbidding the decay. 

Suppose that neutrino masses and mixing arise due to the Majorana mass term of the 
three flavour neutrinos: ^ 

C'^\x) = - - mu'T^iR ui'L + h.c. (36) 

where u^j^ = C^i'il)'^, I = e, ^,t, C being the charge conjugated matrix. We have mw = rriin, 
1,1' = e,fi,T (see, e.g. |13]). The mass term in Eq. fl36l) is diagonalised using the congruent 
transformation: m = U*m'^U\ where m'^ = diag{mi,m2,m^) is a diagonal matrix formed 
by the masses of the Majorana neutrinos Vj and U is the PMNS matrix [^. The effective 
Majorana mass <m> arises in A{PP)o^ from the virtual neutrino propagator (see, e.g. 

!): 

^ = E^S^^ = ^i + ^3 + P5 + ..., (37) 

j 



- rrij 



where 



^l = ^Ef^e>. = ^<"^>, (38) 

^3 = ^Ef^e>,^, etc. (39) 



Here q is the momentum of the virtual neutrino and we have used the fact that <^ 
Typically one has for the average momentum of the virtual neutrino in (/?/?)oj/— decay (see, 
e.g. [70]): ~ (lOMeV)^. As a consequence, the following inequalities hold |P2n+i I ^ 
n = 1,2,.... Usually the terms P3, P5, etc. are neglected in the expression for V. The 
dominant term Pi oc <m>, which leads to A{PP)qi, oc <m>. The factor in Pi gives 
rise to a Coulomb-like potential of interaction between the nucleons exchanging the virtual 
neutrino in the nucleus undergoing (/3/3)oi/— decay. 

Assume now that |<m>| =0. In this case Pi = and the dominant term in the 
expression for V, Eq. fl37|) . will be P3. If |<m>| = is not a consequence of a conservation 
of some lepton charge, we will have P3 7^ and A{PP)qi, 7^ 0, in general. However, unless the 
(/?/?) oj/— decay amplitude receives contributions from mechanisms other than the exchange of 
the light Majorana neutrinos uj, the (/5/3)o,y— decay rate will be extremely strongly suppressed 
due to the fact that [7D] < 10~^^, where we have used nij < 1 eV. Although 

allowed, (/3/3)oy— decay will be practically unobservable if the P3 term in V gives the dominant 
contribution in A{PP)q^. 



^We work in the basis in which the charged lepton mass matrix is diagonal. 
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It is well-known (see e.g. [13]) that |<m>| = |mee|, where m^e is the ee— element of 
the Majorana mass matrix m of neutrinos, Eq. fl36|) . If |<m>| = |mee| = 0, the term 
^^eR ^eL will effectively be "regenerated" at higher orders from the other terms in C^{x), 
Eq. fl36|) . The exchange of virtual Ue mediated by this term will lead to (/5/3)oy— decay. If 
we treat C^^{x) as an "interaction" term[^ and use perturbation theory, the virtual neutrino 
propagator in the (/5/5)ojy— decay amplitude will have, to leading order in the parameters m///, 
the following form: 

P=- — + ..., 40 

gZ gZ 

where 

m = rrief, m*^ m^g + m*^ m^e + rUer m*^ m^e + m^r m*^ m^^ • (41) 

It follows from the expression for the mass parameter rfi that if rriee = 0, we will have m = 
in the following cases [TH [72l [73] : 

i) mef, = rrier = 0, ii) rue^ = m^^ = 0, iii) rrier = m^^, = 0, iv) m^^ = m^^ = m-rr = 0. 
It is easy to see that the four cases in which m = correspond to the conservation of the 
following lepton charges [72]: i) Lf,, ii) — L^, iii) — -^/i, iv) — — L^. In all these 
four cases the (/3/3)oi/— decay is strictly forbidden. However, all four cases are ruled out by 
the existing neutrino oscillation data (see, e.g. [731 El])- Thus, we can conclude that m 7^ 
and therefore A{PP)qi, ^ 0. 

How large can the mass parameter m be? Using the relation m = U* m'^ W and assuming 
that rriee = <m> = 0, it is not difficult to show that 



m 



Y^^l^r (42) 



Thus, we recover the result obtained earlier by expanding the massive Majorana neutrino 
propagators in power series of vnP-jc^: 

1 rh* 

P = - — + ...=P3 + ... (43) 

gZ gZ 

where P3 is given in Eq. (15^ . Therefore the (/?/?)oy— decay will be extremely strongly 
suppressed if = and y4(/3/3)oi/ 7^ is generated at higher order by the Majorana mass 
term, Eq. ([36]) • 



5 Conclusions 

Present and future searches for neutrinoless double beta decay aim at probing lepton num- 
ber violation and the Majorana nature of neutrinos with remarkable precision. A wide 
experimental program is currently under discussion. Experiments with a sensitivity to the 
effective Majorana mass parameter, |<m>|, down to ~ (50 - 10) meV are in a stage of 
preparation or planning and will take place in the future. These experiments will provide 
valuable information on the neutrino masses and the nature of massive neutrinos. 

^In this case vii^(x) should be considered as zero mass fermion fields having the standard zero mass 
fermion propagator. 
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If future (/?/5)o!/-decay experiments show that | < m > | < 0.01 eV, both the IH and the QD 
spectrum will be ruled out for massive Majorana neutrinos. If in addition it is established in 
neutrino oscillation experiments that the neutrino mass spectrum is with inverted ordering, 
i.e. that Am\ < 0, the absence of a signal in neutrinoless double beta decay experiments 
sensitive to |<m>| ~ 10 meV would be a strong indication that the massive neutrinos uj 
are Dirac fermions. At the same time the alternative explanation based on the assumptions 
that the massive neutrinos uj are Majorana particles but there are additional contributions 
to the (/3/3)oy-decay amplitude which interfere destructively with that due to the exchange of 
light massive Majorana neutrinos, would also be possible. However, if Am\ is determined to 
be positive in neutrino oscillation experiments, the upper limit |<m>| < 0.01 eV would be 
perfectly compatible with massive Majorana neutrinos possessing normal hierarchical mass 
spectrum, or mass spectrum with normal ordering but partial hierarchy, and the quest for 
I < m > I would still be open. Under such circumstances the next frontier in the searches for 
(/3/3)oy— decay would most probably correspond to values of |<m>| ~ 0.001 eV. 

Taking |<m>| = 0.001 eV as a reference value, we have investigated in the present 
article the conditions under which |<m>| in the case of neutrino mass spectrum with 
normal ordering would satisfy |<m>| ^ 0.001 eV. We have considered the specific cases of 
i) normal hierarchical neutrino mass spectrum, ii) of relatively small value of the CHOOZ 
angle ^13, as well as iii) the general case of spectrum with normal ordering, partial hierarchy 
and a value of ^^13 close to the existing upper limit. We have derived the ranges of the lightest 
neutrino mass mi and/or of sin^6'i3, for which |<m>| ^ 0.001 eV, and have discussed 
some related phenomenological implications. We took into account the uncertainties in the 
predicted value of |<m>| due to the uncertainties in the measured values of the input 
neutrino oscillation parameters AttIq, Am\ and sin^6'o. For the latter we have used the 
following prospective la errors: 2%, 2% and 4%, respectively. 

In the present analysis we did not include the possible effects of the uncertainty related 
to the imprecise knowledge of the (/3/3)oy— decay nuclear matrix elements. We hope (perhaps 
optimistically) that by the time it will become clear whether the searches for (/3/?)o;y— decay 
will require a sensitivity to values of |<m>| < 0.01 eV, the problem of sufficiently precise 
calculation of the {f3f3)ou— decay nuclear matrix elements will be resolved. 

We have found that in the case of normal hierarchical (NH) neutrino mass spectrum we 
get |<m>| ^ 0.001 eV for sin^6'i3 ^ (0.01 — 0.02) and any value of the relevant Majorana 
phase (difference) 032, provided the currently determined best fit values of the solar and 
atmospheric neutrino oscillation parameters AtUq , Am\ and especially of sin^ 9q, will not 
change considerably in the future high precision measurements (Fig. 3). For < 032 < 7r/2 
one has |<m>| ^ 2.0 x 10~^ eV for any sin^^i3, while if tt/2 < ^32 < 57r/6, we get 
|<m>| ^ 10~^ eV for any sin^ 6^13 allowed at 3a by the existing data. Values of 032 7^ 
in the indicated ranges are required for the generation of the baryon asymmetry of the 
Universe in the "fiavoured" leptogenesis scenario, in which the requisite CP-violation is 
provided exclusively by the Majorana phase (difference) Q!32 [63] . 

We have investigated also the case when sin^ ^13 has a rather small value, sin^ 6*13 ^ 3 x 
10~^, but the neutrino mass spectrum is not hierarchical. We have found that in this case 
one has |<m>| > 10~^ eV for any value of the relevant Majorana phase 021 if the sum of 
neutrino masses satisfies mi + m2 + rris ^ 7 x 10~^ eV. 

In the general case of neutrino mass spectrum with partial hierarchy (i.e. non-negligible 
lightest neutrino mass mi) and sufficiently large sin^^i3, one finds |<m>| > 10^^ eV 
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typically for mi ^ few x 10~^ eV and mi ^ 10~^ eV (Fig. 6). In the second case the 
sum of neutrino masses satisfies mi + 1712 + ^ 7.4 x 10~^ eV. If a future highly sensitive 
(/3/3)oy-decay experiment does not find a positive signal corresponding to |<m>| > 1 meV, 
Majorana neutrinos would still be allowed, but the spectrum would be constrained to be with 
normal ordering and mi to be smaller than ~ 10~^ eV. The CP-parity pattern (+++) will 
be strongly disfavored (if not ruled out) as well. If in addition it is found that sin^ ^13 ^ 0.01, 
mi would be constrained to lie in the interval mi ~ (10-2-10-3) eV (for sin^ ~ 0.32), and 
the CP-parity pattern (+H — ) will also be disfavored. No other future neutrino experiment, 
foreseeable at present, will have the capability of constraining the lightest neutrino mass 
(and the absolute neutrino mass scale) in the meV range. Obviously, the above constraints 
would hold only if massive neutrinos are Majorana particles. If the lightest neutrino has a 
mass in the interval mi ~ (10-^ — 10-^) eV, this can have important effects on the generation 
of the baryon asymmetry of the Universe in the "flavoured" leptogenesis scenario of matter- 
antimatter asymmetry generation [6^ . 

We have provided also an estimate of |<m>| when the three neutrino masses and the 
neutrino mixing originate from neutrino mass term of Majorana type for the (left-handed) 
flavour neutrinos and Yl'j'^^j^ej — 0; but there does not exist a symmetry which forbids the 
{PP)ou-decaj. Our results show that, although in this case the (/?/3)oj/-decay will be allowed, 
the corresponding effective Majorana mass parameter is determined by Y^j ""^^j^ej I 'f ^ where 
q is the momentum of the virtual Majorana neutrino. For the average momentum of the 
virtual neutrino in (/3/3)oi/-decay one typically has (see, e.g. [10]): Ig^l ~ (10 MeV)^. As 
a consequence the contribution to the (/?/3)o,y-decay amplitude y4(/5/?)oi/ due to the light 
Majorana neutrino exchange will be strongly suppressed: |<m>| -C 10-^ eV. Thus, if 
Y^^j ^jU^j = and (/5/3)o!/-decay is observed in an experiment with sensitivity to |<m>| ~ 
10^3 eV, that would imply the existence of contributions to A{PP)qi, due to mechanism(s) 
other than the three light Majorana neutrino exchange. 
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Figure 1: The predicted value of |<m>| as a function of m^^^, obtained using the 2a 
allowed ranges of Am\, Amg, sin^ ^0 and sin^ Oi^. For the NH and QD (and interpolating) 
spectra, the green regions within the black lines of a given type (solid, short-dashed, long- 
dashed, dash-dotted) correspond to the four different sets of CP-conserving values of the two 
phases 021 and a^i, and thus to the four possible combinations of the relative CP parities 
(^2i5'73i) of neutrinos z/12 and z/13: (+1, +1) sohd, (—1,-1) short-dashed, (+1,-1) long- 
dashed, and (—1,-1-1) dash-dotted lines. For the IH spectrum, the blue regions delimited by 
the black solid (dotted) lines correspond to 7721 = +1 {V21 = ~1); independently of rj^i. The 
regions shown in red correspond to violation of CP-symmetry. 
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Figure 2: The predicted value of |<m>| (including a prospective 2a uncertainty) as a 
function of w^^in ^^^^ ~ 0.01. See text for further details. For the NH and QD (and 
interpolating) spectra, the regions within the black lines of a given type (solid, short-dashed, 
long-dashed, dash-dotted) correspond to the four different sets of CP-conserving values of 
the two phases a2i and 031, and thus to the four possible combinations of the relative CP 
parities (7721, "/^si) of neutrinos 1/1^2 and z/i ^: (+1,+1) solid, (—1,-1) short-dashed, (-1-1, —1) 
long-dashed, and (—1, -|-1) dash-dotted lines. For the IH spectrum, the regions delimited by 
the black solid (dotted) lines correspond to 7721 = +1 (??2i = — 1), independently of r^ai. The 
regions shown in red correspond to violation of CP-symmetry. 
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Figure 3: The coloured regions show the values of sin^ ^13 versus Am\ for which | {fTi)\j^j^ < 
1 meV at 1 (2) [3] a (region bounded from below by the dash-dotted (dashed) [dash-double- 
dotted] line for sin^ 9q = 0.32. The error on sin^ 6'i3 is taken to be 0.004 (0.008) in the 
upper (lower) plot. The medium-grey (magenta) region is excluded by the present bound on 
sin^^ia 
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Figure 5: The same as in Fig. 3, but for sin^ 9q = 0.40. 
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Figure 6: The grey (colour online) regions denote the ranges of m^^^ for which |<m>| < 
1 meV and are delimited by thick (thin) lines at 1 (2) a. The CP-conserving patterns are 

indicated by i) solid lines for the case + H — , ii) dashed lines for the H h one, and iii) 

dashed-dotted lines for H . The red triangular region requires CP- violation. The present 

best fit values for Am?, and Ami are used, on 



